Relaxation and flow in linearly sheared two-dimensional foams 
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We probe the relation between shear induced structural relaxation and rheology in experiments 
on sheared two-dimensional foams. The relaxation time tr, which marks the crossover to diffusive 
bubble motion, is found to scale non-trivially with the local strain rate 7 as tr ~ fo ''^7~°'^^. Here 
to is the single-bubble relaxation time, which is four to seven orders of magnitude smaller than 
the inverse strain rate — nevertheless, the flow is not quasistatic. The non-trivial rheology of the 
foam is shown to be intimately linked to the scaling of tr, thus connecting macroscopic flow and 
microscopic bubble motion. 

PACS numbers: 47.57.Bc, 83.60.Fg, 83.80.Iz 

Soft glassy materials, such as foams, emulsions, col- 
loidal glasses and granular media, exhibit highly com- 
plex flows. On the local scale, rearrangements in these 
densely packed systems are erratic due to geometric frus- 
tration. Globally, the relation between strain rate 7 and 
stress T is generally non-linear, often taking a Herschel- 
Bulkley form: r = Ty + cij^ , where Ty denotes the yield 
stress, and where the viscous stress tv = t — Ty scales 
nontrivially with the strain rate 7 

What is the connection between the global rheology 
and the microscopic rearrangements of these materials? 
For non-Brownian systems, the strain rate controls the 
dynamics, therefore a fruitful strategy is to character- 
ize and compare the dependence of both stress and local 
relaxation events on the strain rate. Recent experimen- 
tal and simulation work on sheared foams has uncovered 
rich dynamics on the bubble scale |J, i5|, |6i, i7j, 18] that often 
shows a non-trivial dependence of relaxation time scales 
on the strain rate However, except for ordered 

foams [lo| , a quantitative connection between these time 
scales and the bulk rheology is still lacking. 

In this Letter we establish a direct and quantitative 
connection between the strain rate dependencies of both 
the microscopic relaxation time U, which embodies the 
decorrelation time of the bubble motion, and the macro- 
scopic viscous stress Ty in experiments on a sheared two- 
dimensional, disordered foam. This approach is similar 
in spirit to recent work on glassy systems where this con- 
nection has been probed [ll|, [ij] 

We measure the relaxation dynamics directly by imag- 
ing the bubble motion, while the non-trivial rheology 
of this material is determined from rheometric measure- 
Figure [T] shows examples of the bubble 



FIG. 1: (color online) (a) Schematic top view of the setup, 
with region of interest highlighted and typical flow proflle in- 
dicated by the red curve, (b) Zoom of the highlighted region 
showing short time bubble trajectories. The width of the 
black bar denotes the bin size Ax = 1.6 mm. (c) Example 
of a non-afSne bubble trajectory over a period At = 145 s, 
for 7 = 0.03 s~^. The star and diamond symbol denote the 
starting and endpoint, respectively. 

intermediate times and finally become Gaussian. The oc- 
currence of purely exponential distributions at a sharply 
defined time allows us to extract the relaxation time tr- 
Second, tr is not proportional with the inverse of the 
strain rate, which would be the simplest relation consis- 
tent with dimensional arguments, but instead exhibits 
a non-linear relationship with the inverse strain rate: 
tr = (2.9 ± 0.2) tl'"^-", where u = 0.66 ± 0.05 and to, 
the single bubble relaxation time [l7|, equals 1.3 • 10""* s. 
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trajectories in a top view of our experiment. The aver- 
age fiow profiles exhibit shear bands in our experiment 
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16l | and we employ these to study the statistics 
of bubble displacements over a range of timescales and 
for local strain rates spanning three decades. 

First, we establish that, for a given local strain rate, 
the probability distributions of bubble displacements ex- 
hibit fat tails for short times, develop exponential tails for 



Third, we find that the non-trivial scaling of the vis- 
cous stress with strain rate [1, [3| can be directly re- 
lated to the scaling of the relaxation time: ry = (0.83 ± 
0.05) Gojr, where Go is the static shear modulus and 
the relaxation strain is defined as jr = itr- In addi- 
tion we show that this connection is consistent with a 
non-equilibrium Stokes-Einstein relation. 

Our results firmly link strain rate, relaxation time, and 
stresses in two-dimensional foams. We stress, however, 
that our scenario is not particular to foams, and sug- 
gest to probe similar connections in other sheared, non- 
Brownian systems, such as suspensions, granular media, 
emulsions and microscopic models of these. 

Setup — The measurements are performed in a linear 
shear cell (Fig. la) in which we create a bidisperse 
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FIG. 2; (color online) (a) Evolution of transverse (full symbols) and longitudinal (open symbols) non-afBne displacement 
PDF's at 7 = O.OOSs"^ (■,□) 0.1s; (•,o) 4.1s; (♦,<}) 12.1s. The PDF's are shifted for clarity and normalized by the width 
of the respective distributions. The red curve is a Gaussian distribution. Left inset: The corresponding longitudinal (a) and 
transverse (A) mean square displacements as a function of time. Right inset: the Kurtosis {Asj) / {Asi)^ — 3 for the transverse 
(full symbols) and longitudinal (open symbols) fluctuations. In both insets the three symbols denote the time for which the 
PDF's have been plotted, (b) Mean square displacements (As^) as a function of time at different local shear rates 7 in units 
of s"^ (■) 0.002; (0) 0.009; (a) 0.016; (l>) 0.028; (•) 0.038; (<) 0.11; (T) 0.29. The horizontal line denotes the Lindemann 
criterion (0. 14(d) )^. The inset shows the PDF's of the corresponding mean square displacements at the Lindemann criterion. 



foam layer with bubble diameters of 1.7 and 2.7 mm in 
a soapy solution (5% volume fraction Dawn dishwashing 
liquid, 15% glycerol and 80% demineralized water, vis- 
cosity ?7 = 1.8 • 10"'^ Pa-s, surface tension cr = 28 • 10^'^ 
N/m). The average bubble diameter (d) equals 2.0 ± 0.1 
mm. The bubbles are trapped between the fluid layer and 
a glass plate placed 2.2 mm above the solution. Bubble 
coarsening is negligible over experimental time scales and 
we do not observe rupturing, coalescence or size separa- 
tion of the bubbles. Local shear rates encountered in this 
experiment range from 3 • 10"'* s~* to 0.3 s~^. 

To induce shear, two rotating wheels of 39 cm diameter 
are partially immersed in a soap solution and spaced 7 cm 
= 35 (d) apart. The rotation speed of the wheels controls 
the driving velocity vq of the shearing wall which we vary 
between 0.073 and 2.3 mm/s. The presence of the top 
glass plate causes shear banding [1, Q, and a locally 
varying shear rate 7(0;). We record the bubble motion 
by imaging the monolayer from above at 10 frames per 
second. Custom tracking software allows us to determine 
the bubble positions Si to about 0.005(d). We focus on a 
central region (Fig. la) where average flow transverse to 
the shear direction is absent: {vx{x)) = 0. 

Mean Squared Displacements — In order to relate the 
rearrangement rate of the bubbles and the local strain 
rate, we divide the measurement area into bins of size 
Aa; = 1.6 mm (Fig. lb), and determine for each the lo- 
cal average flow v{x) and shear rate 7 = dv{x)/dx. We 
then determine the non-ajfine bubble tracks Asi(Ai) by 
subtracting the affine mean flow from the bubble trajec- 
tories Si-. Asi(At) = Si{t + A<) - Si{t) - {vi{x))At. An 
example of an erratic bubble trajectory Asi is shown in 
Fig. W[c). We only consider the mean square displace- 



ments on short time scales where y^As|) < (d) , so that 
Taylor dispersion is negligible. 

In Fig. mja), the probability distribution functions 
(PDF's) of both the x and y component of As are shown 
for a given local strain rate 7 at three different times as 
indicated in the left inset of Fig. ^a). The fluctuations 
of the bubble trajectories are isotropic, and we normal- 
ize the widths of all distributions by y^Asf) to high- 
light the qualitative changes in these PDF's with time. 
For short times, they exhibit "fat tails", similar to the 
instantaneous velocity fluctuations observed in a bubble 
raft Then, the PDF's develop exponential tails and 
eventually become Gaussian in the long time, diffusive 
regime. This behavior is also reflected in the evolution of 
the Kurtosis (inset Fig[2Ja)). We note that the time at 
which the PDF's have exponential tails is well defined. 
This behavior should be contrasted to sheared glasses, 
where the PDF's were found to have exponential tails for 
a range of early times before crossing over to Gaussian 
behavior at late times 20j . 
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In Fig. [HJb) we show the longitudinal mean square dis- 
placements for a range of local strain rates as function of 
2l|. We find that the purely exponential PDF's oc- 
(0.14(fi))^. This also holds for the trans- 
"2^ . Interestingly, this corresponds 
to the Lindemann criterion for melting in atomic solids 
and cage breaking in colloidal suspensions ll| . Here the 
relaxation time tr marks the regime where the bubble 
motion crosses over from super-diffusive to diffusive (see 
Fig. [2lb)). Even though tr can thus, in principle, be 
defined from the PDF's or from this crossover, for sim- 
plicity, we define tr as the time where As^ = (0.14(c?))^. 
Note that the non-affine bubble displacements between 
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FIG. 3: (color online) Relaxation time as a function of local 
shear rate for different driving velocities: (■) vq — 0.076 
mm/s; (Q) wo = 0.25 mm/s; (A) vq = 0.62 mm/s; (v) 
Wo = 1.5 mm/s; (♦) vq = 2.3 mm/s. tr is defined as 
{Asl){At = tr) = (0.14(d))^. The line is a fit to the data; 

. on j-0.34 • -0.66 

rearrangements are of similar magnitude (Fig. [T]^c)). 

Relaxation time — We now extract tr for a range of 
local strain rates 7 2]J. Since Brownian fluctuations do 



not play a role, the rearrangements in the foam are en- 
tirely shear induced, and a priori 'y~^ would be a prime 
candidate for setting the relaxation time. In contrast, we 
find that tr = 2.9 tl^^j^", where v w 0.66 ± 0.05, and 
to is the characteristic relaxation time of single bubbles 
r]{d)/a = 1.3 ■ 10~^s — see Fig. |31 This non-trivial scal- 
ing is our central result, and it is far from obvious how it 
can be obtained, if at all, from the local interactions of 
the bubbles. 

Data collapse for mean squared displacements — We 
now use our findings to collapse the mean squared dis- 
placement curves onto a single master curve. Fig. [3Ka) 
illustrates that the system is not quasistatic, since the 
non-affine mean square displacements do not simply scale 
with strain. In fact, Fig.[3]Ja) shows that for a given local 
average strain, Asy is larger in regions where the local 
strain rate is smaller — fluctuations thus increase for 
slower flows. 

Fig. m^b) shows that it is possible to obtain good data 
collapse by rescaling the time axis with the relaxation 
time. We then find a super-diffusive behavior with an 
initial slope « 1.7 for short times. Since large displace- 
ments due to bubble rearrangements are observed even 
on the 0.1 s time scale with which we image the foam, we 
do not see a ballistic regime. At later times the bubbles 
become diffusive and the slope approaches 1. From the 
collapse shown in Fig. SKb) we deduce that the diffusion 
constant Dy scales non-linearly with the shear rate 7 as 
Dy = (0.017 ± 0.002) (d)/ir cx 7". Similar results are 
found for the transverse fluctuations |22] . 



Rheology — We will now establish a surprising relation 
between relaxation times and rheology. Fig. [5] illustrates 
that our 2D foams exhibit a non-trivial rheology, with 
the shear stress t equal to the sum of a yield stress Ty 



and the rate dependent viscous stress Ty = Ci7' 



Re- 



markably, the macroscopic rheological exponent (3 and 
the microscopic relaxation time exponent v are related 
as /9 ~ J/ — 1 . 

This connection between the viscous stress ry and tr 
can be made intuitive by defining a relaxation strain 7^ 
as the product of relaxation time and strain rate. The 
relaxation strain is thus increasing with shear rate, which 
is consistent with observations by Rouyer et al. [23|. To 
convert this relaxation strain into a stress, we multiply 
it with the static shear modulus Go which we measured 
to be 0((jI id)) (see inset Fig. [5]) consistent with Princen 
et al. [2J|. As Fig. [5] shows, Gotrj and the viscous stress 
TV are equal up to a numerical factor of order one: Ty = 
0.83 Gotr'j. Hence the non-trivial global rheology of our 
foams, and in particular the non-trivial power law scaling 
of the viscous stress [sj , is deeply connected to the scaling 
of the relaxation time tr with local strain rate. 

Stokes- Einstein Relation — In previous numerical 
work on sheared foams, a non-equilibrium Stokes- 
Einstein relation was found to relate the product of diffu- 
sion constant, viscosity and bubble size to an energy scale 
obtained from fluctuations in the elastic energy 
We show now that our connection between relaxation 
time and viscous stress may also be viewed from this per- 
spective. Both the diffusion constant Dy and the plastic 
viscosity rip = Ty/7 are strain rate dependent. However, 
the product of these quantities yields an energy scale 
that is constant in the experimentally accessible range 
of strain rates: E — Dyrjp{d) w 7 • 10^'^ cr{d)'^, consistent 
with the simulation results at low strain rates [25i| . Since 
the clastic response of a single bubble can be modelled 
by a spring with a spring constant Aira, the energy scale 
7- 10^"^ cr('^)^ corresponds to compressin g a single bubble 
by « 0.04((i), which appears reasonable [26|. 
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FIG. 4: (color online) (a) The non-affine mean square dis- 
placement (Asy) versus local strain 7 = jAt for different 
local shear rates (same data as Fig. [Sfb)). (b) Same data 
as in (a), but plotted versus dimensionless time At/tr{j) oc 
The dotted line indicates the diffusive limit: 



(Asy) = 2DyAt, with Dy 



(0.017 ±0.002)(d)7fr(7). 
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FIG. 5: (color online) Foam rheology. The stress is plotted 
as a function of shear rate measured in a Couette cell with 
inner radius 5 cm and outer radius 7 cm. The disordered 
foam is floating on the soap solution without the top plate 
and the boundaries have grooves to prevent bubble slip. The 
shear rate at the inner rotating wheel has been determined 
by image velocimetry. (■) total stress. The Herschel-Bulkley 
exponent is P — 0.36 ± 0.05. The curve through the data 
is Go (7y + 0.83 try); (O) viscous stress tv = r — ty where 
TY = 0.22 ± 0.2 Pa. The curve through the data is 0.83 Uj. 
Inset: Oscillatory strain sweep at 0.1 Hz. G' approaches the 
static shear modulus Go at low strain amplitude: Go 7 Pa 
and the static yield strain 7y — 0.031. 

Discussion — Our central result is the scaling of the 
relaxation time scales with strain rate: tr ~ tj~'^7^''. 
Defining quasistatic flows as those for which the precise 
bubble trajectories are independent of the flow rate, the 
nontrivial scaling of tr and Fig. 0^ and [3Jd clearly show 
that our flows are inconsistent with a quasistatic picture, 
even though the global stresses approach rate indepen- 
dence for slow flows. The microscopic scenario which 
emerges is that for increasingly slow flows, the delicate 
balance between elastic and viscous forces [27] causes 
the (relative) bubble motion to become increasingly con- 
centrated in short bursts, such that the contribution of 
viscous forces to the averaged stress vanishes — lead- 
ing to rate independence of the stress [8] . This scenario 
is consistent with the growth of Asy at fixed strain for 
decreasing 7 (Fig. [Hi), and also with the absence of a 
ballistic regime (Fig.lHs): snapshots of the displacement 
fields over a finite time interval reveal strong spatially and 
temporally intermittent behavior, down to the shortest 
timescale that we have probed. 

Finally, it may seem counter-intuitive that our system 
is not quasistatic, given that the single bubble relaxation 
time to [13] is several orders of magnitude smaller than 
the typical inverse strain rate in our system. However, 
the fiuid fiow is concentrated in the thin liquid films sepa- 
rating the bubbles. The shear rate in these films is there- 
fore much closer to 1/to, and the effective viscosity of the 
foam is much higher than the viscosity of the fluid phase 
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